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The distribution of FCoV Types I and II in a Portuguese cat population was studied by a RT- 
PCR assay targeting the 3'-end of the viral RNA. For a period of 3 years, 120 samples were 
collected and 57 were found positive for FCoV RNA. Within the positive samples the 
presence of FCoV Type I was found in 79%. Type II was only detected in 3.5% in animals 
with Feline Infectious Peritonitis. The remaining 17.5% could not be differentiated. These 
viral sequences, comprising a region within gene S were further subjected to a 
heteroduplex mobility assay (HMA) detecting the presence of viral quasispecies in 17% of 
the samples. Phylogenetic analysis for FCoV Type I revealed high genetic diversity 
between the Portuguese sequences and other previously characterized strains, while 
Type II tree showed a higher genetic homogeneity. This study confirmed the presence of 
FCoV Types I and II circulating in Portugal and detected high genetic diversity between 
circulating strains suggesting that the virus persists within the host as mixed viral 
populations. 

© 2009 Elsevier B.V. All rights reserved. 


1. Introduction 

Feline Coronaviruses (FCoVs) are members of genus 
Coronavirus, Family Coronaviridae. Two variants of FCoV 
have been described, Feline Enteric Coronavirus responsible 
for an enteric disease and Feline Infectious Peritonitis 
Virus (FIPV) the etiological agent of FIP, a fatal disease of 
domestic and wild cats. FIP development in a small 
percentage of animals infected with FCoV is related to the 
occurrence of genomic mutations supported by the close 
genetic relationship between the higher virulent virus and 
the parental FCoV (Herrewegh et al., 1997; Vennema et al., 
1998; Vennema, 1999). Additionally, FCoVs were classified 
as Types I and II according to their serological and genomic 
features. Type I is strictly feline. Type II FCoV, although 
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sharing with Type I a high degree of similarity in the 3'-end 
of the genome, is antigenic and genetically more related to 
Canine Coronavirus (CCoV) (Hohdatsu et al., 1991), 
suggesting that Type II resulted from a homologous and 
independent recombination event within the S and the M 
genes of Type I FCoV and CCoV (Herrewegh et al., 1998). 
The current information on the geographical distribution 
of FCoV indicates a high incidence of Type I in Europe, 
Japan and USA (Addie et al., 2003; Hohdatsu et al., 1992). 
However in Japan and Taiwan detection of FIP was 
associated with Type II FCoV (Hohdatsu et al., 1992; Lin 
et al., 2009). Presently there is no available information 
about the molecular epidemiology of FCoV in Portugal. Due 
to the pathogenic impact of this virus in the feline health it 
was our purpose to investigate the presence of FCoV Types 
in the Portuguese cat population using a RT-PCR-nested 
PCR assay previously reported (Addie et al., 2003) and to 
assess its genetic diversity by heteroduplex mobility assay 
(HMA) and phylogenetic analysis. 
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2. Materials and methods 

2.2. Sampling collection 

The biological materials included rectal swabs, 
stool, plasma, thoracic and ascite fluids collected from 
animals attending the Veterinary Hospital of the Faculty 
of Veterinary Medicine, Lisbon. 

2.2. Detection of FCoV RNA 

FCoV RNA was extracted using the QIAmp Viral RNA or 
Qiagen RNeasy mini kits (Qiagen), following the manu¬ 
facturer’s instructions. A fragment of 177 bp located in the 
3'-UTR was amplified by RT-PCR followed by a nested PCR 
using two primer pairs described by Herrewegh et al. 
(1995). The 25 pX reactions were amplified with one step 
RT-PCR Fidelity Taq 2x Master Mix (GE Healthsciences), 
50pmol of each primer and 10 pX of RNA sample. The 
nPCR was performed with PCR Fidelity Taq 2x Master Mix 
(GE Healhsciences), 50 pmol of each primer and 5 |ulL of the 
previous amplification in a 25 |mL reaction 

2.3. Differentiation between Types I and II FCoV 

Positive FCoV RNA samples were subjected to a RT-PCR 
assay to amplify a 376-283 bp fragment located in the 3'- 
region of the spike gene (S), followed by a nPCR to amplify a 
360-218 bp internal fragment, for differentiation of Type I 
from Type II FCoV according to the fragment length (Addie 
et al., 2003). The amplification was performed with one 
step RT-PCR Fidelity Taq 2x Master Mix (GE Healths¬ 
ciences), 5 |ulL of RNA sample and primers Iffs, Icfs, Iubs in a 
25 pX reaction, followed by a nested PCR with PCR Fidelity 
Taq 2x Master Mix (GE Healthsciences), primers nlffles, 
nlcfs, nlubs (Addie et al., 2003) and 2.5 pX of the previous 
reaction in 25 pX final volume. The amplicons were 
analyzed in a 2.5% agarose gel stained with 0.5 p,g/mL of 
ethidium bromide in 1 x TBE. 

2.4. Heteroduplex mobility assay (HMA) and cloning 

The amplicons were subjected to a HMA before cloning 
(Delwart et al., 1993). Briefly, 10 |uX of the second round S 
amplicons were combined with 1 pX of lOx annealing 
buffer (1 M NaCl, 100 mM Tris pH 7.8, 20 mM EDTA), 
heated to 94 °C/5 min, cooled on ice, added of 1 |uX of 
loading buffer (50% glycerol, 0.01% of xyleno cyanol) and 
loaded in a 5% polyacrylamide gel (29.2:0.8 acrylamide:- 
bisacrylamide) poured into plates of 15 cm x 15 cm. The 
electrophoresis was carried at a constant 150 V/3 h in a 
vertical gel apparatus (Biorad) in lx TBE. After electro¬ 
phoresis the gel was stained with 0.5 p.g/mL of ethidium 
bromide in 1 x TBE and observed in a UV transiluminator. 
The visualized homoduplexes bands were cuted, eluted at 
56 °C/5 min in 200 |uX H 2 0 of which 1 pT was re-amplified 
using the nPCR primers set described above. The resulting 
amplicons were cloned in pGEM®-T Easy Vector System 
(Promega) according to the manufacture’s conditions and 
sequenced by STAB Vida Lda, Portugal. The nucleotide 
sequences were submitted to Genbank and were assigned 
the accession numbers EU327694-EU327752. Sample 


identification was included in Table 1. FCoV. Types I and 
II assignment based on the amplicons length were in 
accordance with the nucleotide sequence of each sample. 

2.5. Phylogenetic analysis 

The nucleotide sequences were aligned (PILEUP from 
GCG Wisconsin Package v. 10.0), including sequences of 
FCoV Type I (KU_2 AN-D32044; UCD1 AN-AB088222; Black 
AN-AB088223; C2_0301UK AN-AY159737; J8_0196UK AN- 
AY159758; FQ1_1099UI< AN-AY159754; H9_l 100UK AN- 
AY159770; FUL1191UK AN-AY159735) and FCoV Type II 
(79_1683 AN-X80799; 79_1146 AN-X06170; G6_1095 AN- 
AY159784). Due to its recombinant status FCoV Type II 
sequences were aligned and submitted to phylogenetic 
analysis separately from Type I sequences. 

The phylogenetic relationship between sequences was 
inferred by genetic distances between pairs of sequences 
using Kimura’s two parameter method (DNADIST, Phylip 
Package). The values of transition/transversion ratio and 
nucleotide frequencies were calculated from the data sets 
(TREE-PUZZLE) (Strimmer and von Haeseler, 1997). The 
resulting trees were constructed by Neighbor Joining 
(NEIGHBOR, Phylip Package) and the branching order 
reliability was evaluated by bootstrap analysis of 1000 
replicates (SEQBOOT, Phylip Package). To maximize 
the phylogenetic information of the nucleotide align¬ 
ment closely related sequences were excluded but were 
indicated between brackets following the sequence 
chosen to represent each set of viruses as summarized; 
65LVL05 [74LVL05, 81LVL05, 86LVL05, 87LVL05]; 
14LVI_05 [ 10LVI_05a]; 19LVI_04a [ 19LVI_04b]; 
58LVL05 [59LVI_05a/b]; lLg_06 [lFg_06a]. 

3. Results 

3.2. Detection of Types I and II FCoV 

For 3 years 120 samples were collected from domestic 
cats to detect FCoV RNA due to clinical suspicion of FIP and 
to confirm viral excretion from healthy animals. Fifty- 
seven samples were found FCoV positive (47.5%) and 
afterwards processed for Types I and II differentiation. Out 
of those 57 positive samples, 45 (79%) were classified as 
Type I FCoV and were obtained from healthy animals or 
with FIP. Only two samples (3.5%) from animals with FIP, 
were classified as Type II. No amplification was obtained 
from ten samples (17.5%). 

3.2. Detection of viral quasispecies 

The presence of viral quasispecies within the animal, 
was evaluated by HMA for each sample (Fig. 1). This assay 
detects the presence of fully (homoduplexes) and partially 
(heteroduplexes) complementary DNA chains formed 
between highly related sequences during the process of 
denaturing and annealing. Homoduplexes molecules tend 
to migrate faster, while heteroduplexes show a reduced 
mobility in a non-denaturing PAGE. The observation of two 
or more homoduplexes bands, plus the corresponding 
heteroduplexes bands indicates the presence of different 
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Table 1 

Identification and origin of the Portuguese FeCoV sequences, amplified from healthy and FIP (S) animals. The genetic diversity of the amplicons was 


evaluated by heteroduplex mobility assay and each positive sample was 

re-amplified. The number of corresponding sequences is indicated (a-c). 

Identification 

Heteroduplexes 

Origin 


2004 




13LVL2004 

Positive (a;b) 


Lisbon, Portugal 

19LVL2004 

Positive (a;b) 


Lisbon, Portugal 

20LVL2004 

Negative 


Lisbon, Portugal 

22LVL2004 

Positive 

Palmela, Portugal 

Sharing household 

23LVL2004 

Negative 

Palmela, Portugal 


25LVL2004 

Positive (a;b) 

Palmela, Portugal 


2005 




3LVL2005 

Negative 


Lisbon, Portugal 

10LVL2005 

Positive (a;b) 

Lisbon, Portugal 

Sharing household 

11LVL2005 

Positive (a;b) 

Lisbon, Portugal 


14LVL2005 

Negative 

Lisbon, Portugal 


27LVL2005 

Negative 


Lisbon, Portugal 

33LVL2005 

Positive (a;b) 

Rio Mouro, Portugal 

Multi-cat environment 

37LVL2005 

Negative 


Lisbon, Portugal 

38LVL2005 

Negative 


Lisbon, Portugal 

48LVL2005 

Negative 


Lisbon, Portugal 

50LVL2005 

Negative 


Lisbon, Portugal 

58LVL2005 

Negative 


Lisbon, Portugal 

59LVL2005 

Negative 


Lisbon, Portugal 

71LVL2005 

Negative 


Lisbon, Portugal 

80LVL2005 

Negative 


Vila Nova de Gaia, Portugal 

65LVL2005 (S) 

Negative 

Lisbon, Portugal 

Sharing household 

74LVL2005 

Negative 

Lisbon, Portugal 


81LVL2005 

Negative 

Lisbon, Portugal 


86LVL2005 

Negative 

Lisbon, Portugal 


87LVL2005 

Negative 

Lisbon, Portugal 


88LVL2005 

Negative 

Lisbon, Portugal 


90LVL2005 

Negative 

Lisbon, Portugal 


93LVL2005 

Negative 

Lisbon, Portugal 


94LVL2005 

Negative 

Lisbon, Portugal 


96LVL2005 

Negative 

Lisbon, Portugal 


99LVL2005 

Positive (a;b) 


Lisbon, Portugal 

102LVL2005 

Negative 


Lisbon, Portugal 

2006 




12LVL2006 

Negative 


Lisbon, Portugal 

13LVL2006 

Negative 


Lisbon, Portugal 

15LVL2006 

Negative 


Lisbon, Portugal 

19LVL2006 

Negative 


Lisbon, Portugal 

20LVL2006L 

Negative 


Lisbon, Portugal 

20LVL2006PBMC 

Negative 


Lisbon, Portugal 

20LVL2006Z 

Negative 


Lisbon, Portugal 

24LVL2006 

Negative 


Lisbon, Portugal 

24LVL2006D1 

Negative 


Lisbon, Portugal 

1FG_2006 (S) 

Positive (a;b) 


Lisbon, Portugal 

1LG_2006 

Negative 


Lisbon, Portugal 

28LVL2006 

Negative 


Lisbon, Portugal 

29LVL2006 

Negative 


Lisbon, Portugal 

32LVL2006 (S) 

Positive (a;b,c) 


Lisbon, Portugal 

38LVL2006 

Negative 


Lisbon, Portugal 


DNA populations within the sample, named viral quasis¬ 
pecies and mirror the in vivo genetic diversity of the virus 
(Delwart et al M 1993). 

Within the 47 S gene amplicons 37 samples (78.7%) 
showed no genetic diversity whereas in 10 samples (21.3%) 
two homoduplexes bands were observed and individually 
amplified. Comparison of the nucleotide sequences of the 
heteroduplexes positive samples revealed similarities of 
90-99% along the 320nt in Type I and 170nt in Type II FCoV. 
Although very similar, six sequences (13LVI_04b; 
19LVI_04b; 10LVI_05a; llLV_05a; 33LVI_05b; 99LVI_05b) 
showed a 5'-end variation length of 28-56 nucleotides 
justifying the homo/heteroduplexes formation. Sequences 


33LVI_05a/b, 25LVI_04a/b, 10LVI_05a/b and llLVI_05a/b 
were amplified from animals with multi-cat environment 
(Table 2). 

3.3. Phylogenetic analysis 

The genetic relationships between the Portuguese 
sequences of FCoV and previously assigned sequences 
from isolates of Type revealed a heterogeneous cluster¬ 
ing of sequences with several tree branches supported 
by high bootstrap values (Fig. 2). Sequences from 
distinct viral quasispecies obtained from the same 
animal as well as sequences from cats sharing the same 





166 


A Duarte et al. / Veterinary Microbiology 138 (2009) 163-168 


123456789 10 


* 

* 


* * 

* * 


¥4 



Swuu 



+ + + + + + 


Fig. 1. Heteroduplex mobility assay of S amplicons from samples 
10LVL05 (1); 11LVL05 (2); 3LVL05 (3); 25LVL04 (4); 23LVL04 (5); 
22LVL04 ( 6 ); 20LVL04 (7); 10LVL04 ( 8 ); 13LVL04 (9); 12LVL04 (10), 
included in Type I FCoV. Lanes 1,2,4, 6 , 8 and 9 (+) show the presence of 
two homoduplexes bands (-»), corresponding to fully complementary 
DNA chains. The heteroduplexes bands (*) consistent to partially 
complementary DNA chains with reduced mobility are clearly visible 
in lanes 2, 6 , 8 and 9. In lanes 3, 5 and 7 no homo/heteroduplexes are 
present. The photograph of the ethidium bromide-stained gels is 
displayed as reversed image. 


household were included in branch A indicating the 
presence of different viral strains within each confined 
cat population and within the same individual host. 
Sequences 83 LVL05, 88LVL05, 90LVL05, 94LVL05 and 
96LVL05 came from cats sharing the same household 
than the source animal of sequence 65LVL05 [74LVL05, 
81LVL05, 86LVL05, 87LVL05]. However sequence 
65LVL05 was positioned outside group A, in a separate 
branch with 80% of bootstrap value. The above-men¬ 
tioned sequences were obtained from a multi-cat 
household where a clinical FIP episode (65LVL05) had 
been reported. Although all the other animals were FCoV 
positive, none developed FIP. Sequence llLVI_05b was 
also included in branch A but sequence llLVI_05a, 
obtained from the same cat was positioned in branch B. 
Sequence 25LVI_04a was grouped in branch A though it 
was obtained from the same animal as sequence 
25LVI_04b placed in a different branch. These two viral 
quasispecies (25LVI_04a/b) came from the same house¬ 
hold of the source cats of sequences 22LVI_04a/b and 
23LVL04, which in this analysis were found distributed 
in outlier branches supported by 100% of bootstrap 
value. Sequences 20LVL06Z (stool), 20LVL06PBMC 
(blood) and 20LVL06L (ascites) amplified from the same 
animal, were placed in two separate branches, one with 
100% bootstrap value (20LVL06PBMC and 20LVL06L). 

Type II FCoV phylogenetic tree included previously 
assigned FCoV Type II sequences and two Portuguese 
sequences obtained during this study, one positive to viral 
quasispecies (32LVL06; 32LVI_06a/b). All the sequences 
clustered in one branch, although FIPV 79-1146 and 79- 
1683 grouped together with 100% of bootstrap value. Due 
to the low number of nucleotides considered in the align¬ 
ment (170b), a likelihood mapping analysis was preformed 
(TREE-PUZZLE) to assess its phylogenetic content, revealing 
62.86% of fully resolved, 0% of partially resolved and 37.14% 
of unresolved quartets demonstrating therefore its reliable 
use to estimate phylogenetic relationships. 


Table 2 

Genetic diversity between amplicons subjected to HMA. The length and 
the percentages of genetic similarity between homoduplexes are shown. 


Identification 

Sequence 

Similarity (%) 

FCoV 

13LVI_04a 

320 bp 

317/320 (99%) 

Type I 

13LVI_04b 

348 bp 

Extra 5' 28nt 


19LVI_04a 

320 bp 

320/320 (100%) 


19LVI_04b 

348 bp 

Extra 5' 28nt 


22LVI_04a 

22LVI_04b 

347 bp 

347 bp 

344/347 ( 99 %) 


25LVI_04a 

25LVI_04b 

320 bp 

320 bp 

290/320 (90%) 


10LVI_05a 

348 bp 

316/320 (98%) 


10LVL05b 

320 bp 

Extra 5' 28nt 


11 LVI_05a 

348 bp 

303/320 (94%) 


llLVL05b 

320 bp 

Extra 5' 28nt 


33LVI_05a 

320 bp 

319/320 (99%) 


33LVI_05b 

348 bp 

Extra 5' 28nt 


99LVI_05a 

320 bp 

319/320 (99%) 


99LVI_05b 

376 bp 

Extra 5' 56nt 


1 FLVI_06a 

1 FLVI_06b 

320 bp 

320 bp 

299/320 (93%) 


32LVI_06a 

32LVI_06b 

170 bp 

170 bp 

168/170 (98%) 

Type II 


4. Discussion 

Until now, no previous study on the presence and 
prevalence of FCoV Types I and II in domestic cats had 
been done in Portugal. In our analysis 47.5% of the 
samples were positive in the RT-PCR assay (57/120) and 
were obtained either from healthy animals or from cats 
with FIP. Moreover, the vast majority of the amplicons 
was classified as Type I. Type II was found only in two 
samples both from cats with FIP. The values found for 
FCoV Types I and II were in accordance with the 
distribution reported in Europe and USA (Addie et al., 
2003; Benetka et al., 2004) although in Japan and Taiwan 
a higher prevalence of FCoV Type II was observed in 
animals with FIP (Hohdatsu et al., 1992; Lin et al., 2009). 
In our study even though the two FCoV Type II sequences 
were obtained from sick cats, the low incidence of this 
viral biotype does not allow any relationship between 
FCoV II and symptomatic disease. 

The evaluation of FCoV genetic diversity revealed 21% 
of the amplicons positive for viral quasispecies. These 
viral variants were found within the same animal but also 
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Fig. 2. Phylogenetic analysis of FCoV Types I and II sequences based on the 3'-end of the S gene. The trees were constructed by Neighbor Joining based on the 
genetic distances between the nucleotide sequences, calculated by Kimura two-parameter using a transition/transversion rate of 3.23 for Type I sequences 
and 9.99 for Type II sequences. A boostrap analysis of 1000 replicates was performed and values above 80% are indicated on the branches. Branches A-C 
include FCoV sequences with the lowest value of nucleotide substitutions within each branch. Sequences marked with an asterisk (*) were amplified from 
animals with FIP. 


in samples collected from animals sharing a multi-cat 
environment. This finding supports the viral internal 
mutation theory (Dye and Siddell, 2007) and confirms the 
animal exposure to divergent viral strains. The signifi¬ 
cance of FCoV quasispecies had already been studied 
(Battilani et al., 2003) suggesting a correlation between 
high viral genetic diversity within the host and develop¬ 
ment of disease. However from our data it was not 
possible to correlate viral diversity with virus subtype or 
disease therefore a more extensive study should be 
performed to assess this connection. 

The phylogenetic analysis of FCoV Types I and II 
sequences revealed an assorted distribution of sequences 
in Type I tree, showing high genetic diversity of the FCoV 
Portuguese sequences, its divergence toward European 
and Japanese sequences (Fql_1099UK, H9_1100UK, and 
strains KU2, Black and UCD1) confirming the previous 
findings. 

Interestingly, FCoV Type II phylogenetic tree although 
with a limited number of Portuguese and previously 
characterized Type II sequences, showed a homogeneous 


distribution of sequences as also reported by Lin et al. 
(2009) in a similar study, indicative of high similarity 
between distantly related sequences (UK; Netherlands; 
Portugal). The genetic diversity of the Portuguese FCoV 
sequences confirmed and strengthens the evidence that 
this virus persists within the host as mingled viral 
populations whose interaction with the animal may or 
not result in severe disease. 

In conclusion this study provided information on the 
FCoV molecular epidemiology in Southern Europe which 
was unavailable until now and revealed high genetic 
diversity between the Portuguese and previously char¬ 
acterized FCoV sequences. 
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